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Abstract 


The  adsorption  of  dimethyl  methylphosphonate  (DMMP)  on  the  (010)  surfaee  of  anatase 
TiOi,  whieh  is  isostructural  with  the  (100),  has  been  studied  using  density  funetional  theory  and 
two-dimensionally-periodie  slab  models.  The  experimentally-observed  faeeting  of  this  surfaee 
has,  for  the  first  time,  been  ineluded  in  the  modeling.  The  relaxations  of  bare  surfaees  both  with 
and  without  faeeting  are  similar,  leading  to  an  atomie-seale  roughening  due  to  inward  (outward) 
displaeement  of  five-fold-eoordinated  Tisc  (six-fold-coordinated  Tec)  atoms  together  with  out¬ 
ward  displacement  of  three-fold-coordinated  Osc  atoms.  Molecular  adsorption  occurs  by  forma¬ 
tion  of  a  Tisc — 0=P  dative  bond  with  one  or  more  C-H — 02c  bonds  between  CH3  groups  and  un¬ 
saturated,  two-fold-coordinated  (02c)  sites.  The  energies  for  molecular  adsorption,  obtained  us¬ 
ing  the  B3LYP  functional,  are  virtually  identical  (about  -21.0  kcal/mol)  for  the  two  surfaces  and 
are  also  close  to  those  found  elsewhere  for  the  rutile  (110)  and  anatase  (101)  surfaces.  A  possible 
first  step  in  the  dissociative  adsorption  of  DMMP  has  also  been  modeled  and  is  found  to  be  ther¬ 
modynamically  favored  over  molecular  adsorption  to  a  degree  which  depends  on  faceting. 
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1.  Introduction 


The  surface  structure  and  properties  of  Ti02,  in  both  the  rutile  and  anatase  forms,  have  been 
intensively  studied  both  experimentally  and  computationally  [1],  For  either  form,  the  adsorption 
of  a  wide  variety  of  molecules  has  been  investigated  on  several  different  crystallographic  planes. 
An  especially  important  aspect  is  adsorption  on  the  isostructural  (100)  and  (010)  surfaces  of  ana¬ 
tase'.  Many  of  the  technological  applications  for  Ti02,  in  which  the  material  is  used  as  a  catalyst, 
a  photocatalyst  [2]  or  a  catalyst  support,  involve  high-surface-area  (HSA)  powders  prepared  by 
burning  TiCU  in  O2.  The  Degussa  P25  powder,  which  has  been  widely  used  in  surface-chemistry 
experiments  on  Ti02,  is  an  example  of  such  material.  Typically  these  powders  are  ~80%  anatase, 
and  experiment  [3-5]  shows  that  in  this  case  the  (010)  is  the  dominant  surface.  The  (100)  is  also 
known  to  make  a  major  contribution  to  the  surface  area  of  anatase  nanocrystals  [6]  and  to  be  es¬ 
pecially  active  in  photocatalysis  [7].  Recent  work  by  Dzwigaj  et  al.  [8]  has  clearly  shown  that  the 
morphology  of  anatase  powder  (i.e.,  the  distribution  of  exposed  crystal  planes)  depends  on  the 
method  of  fabrication  and  treatment.  However,  in  general  the  (100)  surface  appears  to  have  a  sig¬ 
nificant,  if  not  always  dominant,  presence  in  many  forms  of  anatase  powder. 

Experimental  data  for  this  surface  are  scarce,  due  to  the  difficulty  in  obtaining  anatase  sam¬ 
ples  suitable  for  single-crystal  studies.  However,  several  computational  studies  of  adsorption 
have  been  reported,  focusing  on  atomic  H  [9-11],  NH3  and  NO  [12],  CO2  [13],  H2O  [12,14], 

O2  [14],  chlorobenzene  [15],  chlorophenol  [16],  N  heterocycles  [17],  acetic  acid  [18]  and  oxalic 
acid  [19].  In  most  cases  unsaturated,  five-fold-coordinated  (Tisc)  sites  are  active  as  Lewis-acid 
adsorption  sites.  In  particular  the  adsorption  behavior  of  the  bidentate  reagent  oxalic  acid  on  ana- 

'  The  designations  (100)  and  (010)  will  be  used  interchangeably  since  the  two  anatase  surfaces 
are  isostructural,  as  are  the  (101)  and  (011).  Some  of  the  literature  on  anatase  powder  refers  to 
the  (010),  and  that  labeling  is  used  here  except,  in  the  Introduction,  when  citing  works  which  re¬ 
fer  explicitly  to  (100). 
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tase  (100)  differs  from  that  on  rutile  (110)  due  the  differing  loeal  arrangement  of  Tisc  surfaee 
sites  [19].  Unsaturated,  two-fold-eoordinated  (02c)  sites  ean  also  be  involved  in  hydrogen  bond 
(H-bond)  formation.  The  effects,  on  the  structure  of  the  (100)  and  other  anatase  surfaces,  of  hy- 
droxylation  and  of  steady-state  ambient  pressures  of  Hi,  HiO  and  HiS  have  been  thoroughly  in¬ 
vestigated  in  a  series  of  density  functional  theory  (DFT)  studies  by  Arrouvel  et  al.  [20,21].  These 
studies  have  provided  important  insight  into  the  properties  of  real  TiOi  catalysts. 

Figure  1  shows  schematic  models  for  the  (010)  surface  which,  except  for  the  (101),  has  the 
lowest  surface  energy  (g)  for  anatase  [20,22].  Figure  la  shows  the  ideally-terminated  (unrelaxed) 
surface,  and  Fig.  lb  shows  the  (1x2)  faceted  surface  (also  unrelaxed)  observed  [23]  in  scanning 
tunneling  microscopy.  Faceting  of  the  (100)  surface  exposes  the  more  stable  (101),  which  is  iso- 
structural  with  the  (Oil).  Computation  [20,22]  shows  g  =  0.44  and  0.53  J/m^  for  the  relaxed 
(101)  and  (100)  surfaces,  respectively,  vs.  0.31  J/m  for  rutile  (110).  The  effect  of  faceting  on  the 
anatase  (100)  g  is  unknown  at  present.  All  (100)  surface  Ti  sites  are  Tisc  while  half  the  surface  O 
atoms  are  Oic  and  half  are  saturated,  three-fold-coordinated  Osc  sites.  With  equal  numbers  of  sin- 
gly-unsaturated  Tisc  and  02c  sites  the  surface  is  auto-compensated,  meaning  that  the  excess  elec¬ 
tron  density  on  a  Tisc  is  transferred  to  an  02c  to  maintain  a  closed-shell  configuration  on  both. 
Hence,  there  is  no  indication  [10,17]  of  surface  states  in  the  band  gap  of  the  relaxed  but  unfacet¬ 
ed  surface.  Adsorption  on  the  faceted  surface,  which  has  not  to  our  knowledge  been  studied  pre¬ 
viously,  should  be  especially  interesting.  For  an  appropriately-chosen  molecule  the  structure  of 
this  surface,  with  rows  of  Tisc  sites  set  below  and  in  close  proximity  to  rows  of  02c  sites,  holds 
out  the  possibility  of  both  dative  bonding  to  a  Tisc  and  H-bonding  to  multiple  02c  sites. 

The  purpose  of  this  work  is  to  examine  the  adsorption  of  a  fairly  complex  molecule, 
(CH30)2(CH3)P=0  (dimethyl  methylphosphonate,  DMMP),  on  the  anatase  (010)  surface  both 


4 


with  and  without  faceting.  This  is,  to  our  knowledge,  the  first  computational  study  of  the  faeeted 
anatase  (010)  surfaee  with  or  without  adsorbates.  Our  interest  in  DMMP  stems  from  the  use  of 
this  speeies  as  a  simulant  for  chemical  warfare  agents  (CWAs)  sueh  as  Sarin  and  the  demonstrat¬ 
ed  utility  of  Ti02  in  the  eatalytie  destruetions  of  CWAs  (refs.  24-28  and  works  eited).  The  grow¬ 
ing  interest  in  HSA  Ti02  powders  for  this  applieation,  and  the  wide-spread  use  of  DMMP  in  ex¬ 
perimental  work  relating  to  CWAs,  makes  sueh  a  study  timely.  For  two  reasons,  the  present  work 
foeuses  on  adsorption  of  DMMP  on  bare  surfaees,  free  of  defeets  or  eo-adsorbed  speeies.  Firstly, 
the  adsorption  of  DMMP  on  Ti02,  even  under  these  ideal  eonditions,  remains  eontroversial.  See- 
ondly,  the  experimental  data  for  this  system  [29-32]  with  whieh  we  will  eompare  our  results 
were  all  obtained  in  vacuum  or  in  ultra-high  vacuum  (UHV). 

2.  Computational  and  modeling  details 

2.1  Geometry  optimization  and  energy  calculations 

All  calculations  were  done  using  two-dimensionally-periodie  slab  (2-DPS)  models  and  the 
CRYSTAL  06  eode  [33-35]  whieh  employs  Bloeh  funetions  construeted  from  loealized  Gaussian 
basis  sets.  To  reduee  the  eomputational  eost,  geometry  optimization  was  done  at  the  restrieted 
Hartree  Foek  (RHF)  level.  This  has  previously  been  shown  [36,37]  to  give  results  for  the  bulk 
and  surfaee  struetures  of  Ti02  whieh  are  very  elose  to  those  obtained  using  DFT.  Optimization  at 
the  RHF  level  has  also  been  shown  to  give  reliable  results  for  moleeular  geometries  [38].  For  Ti, 
the  Hay-Wadt  small-eore  effeetive-eore  pseudopotential  (ECP)  was  used  with  the  41 1(31  ld)G 
basis  set  previously  optimized  for  perovskite  titanates  [33,39].  For  O,  a  6-31G(d)  all-eleetron  ba- 
sis  set  [40]  was  used  with  the  exponents  of  the  outermost  shells  re-optimized  for  bulk  rutile.  The 
optimized  parameters  were  0.27281  and  0.31201  Bohr'  ,  respeetively,  for  the  sp  and  d  shells.  The 

^  Basis-set  optimization  was  done  using  CRYSTAL  03  with  the  LoptCG  script  (C.  Zicovich- 
Wilson;  Universidad  Autonoma  del  Estado  de  Morelos;  Cuernavaca,  Mexico). 
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bulk  lattice  constants  and  atom  positions  optimized  using  this  procedure  show  good  agreement 
with  experiment  (Table  1).  In  subsequent  geometry  optimizations,  the  lattice  constants  were  fix¬ 
ed  at  these  values  while  the  positions  of  all  Ti  and  O  atoms  (as  well  as  those  of  the  DMMP  at¬ 
oms)  were  allowed  to  vary.  For  DMMP,  6-31G(d,p)  basis  sets  from  a  standard  database  [40] 
were  used  without  modification. 

After  geometry  optimization,  a  single-point  calculation  of  the  total  energy  was  done  using 
DFT  with  all-electron  basis  sets  and  the  B3LYP  functional,  which  has  been  shown  [37,41]  to 
give  reliable  results  for  the  bare  surfaces  of  both  rutile  and  anatase.  A  discussion  of  the  use  of 
B3LYP  in  describing  the  adsorption  of  DMMP  is  given  elsewhere  [42].  For  Ti02  the  basis  sets 
are  designated  86-41  l(d31)G  for  Ti  and  8-41  IdlG  for  O,  while  the  6-31G(d,p)  basis  sets  were 
again  used  for  DMMP.  For  Ti  and  O  the  outermost  sp  and  d  shells  have  been  optimized  [41]  for 
TiOi,  and  the  complete  basis  sets  are  obtained  by  combining  these  with  the  inner  shells  given 
elsewhere  [33].  Preliminary  calculations  using  somewhat  larger  6-31  lG(d,p)  basis  sets  [40]  for 
DMMP  led  to  difficulties  in  converging  the  self-consistent  field  (SCF).  This  may  result  from  the 
small  exponent  in  the  outermost  sp  shell  of  the  P  atom  (0.0684930  Bohr'  )  which  can  cause  lin¬ 
ear-dependency  problems  in  calculations  for  periodic  structures  [33]. 

The  bulk  calculations  used  an  (8x8x4)  sampling  grid,  giving  a  k-point  spacing  of  ~0.1  atom¬ 
ic  units  (au).  Calculations  for  the  (4x1)  and  (4x2)  supercells  used  in  adsorption  (see  Sec.  2.2) 
were  done  with  (1x2)  and  (1x1)  grids  respectively,  giving  a  k-point  spacing  of  ~0.2  au.  Reducing 
the  k-point  spacing  to  0.1  au  for  the  (4x1)  supercell  had  a  negligible  effect  on  the  energy.  Fol¬ 
lowing  Labat  et  al.  [37],  an  "extra-large"  (75,974)  integration  grid  was  used  for  all  DFT  calcula¬ 
tions.  In  CRYSTAL,  truncation  of  the  sums  of  Coulomb  and  exchange  terms  in  the  Fock  matrix  is 
determined  by  five  overlap  criteria  [33,35]  (T1-T5).  These  were  set  at  10'^  for  T1-T4  and  10"'^ 
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for  T5.  A  high  numerical  accuracy  in  integration  is  generally  required  for  Ti02  [37,41]  and  is 
further  necessitated  by  the  small  exponent  (0.09983  Bohr'  )  in  the  outermost  5p-shell  of  the  P 
6-31G(d)  basis  set.  In  the  multipolar  expansion  zone  [33-35]  a  maximum  order  of  L=  6  was  used. 
The  geometry  convergence  criteria  were  1x10'  Hartree  for  energy,  1.8x10'  Bohr  for  maximum 
displacement  and  4.5x10'"^  Hartree/Bohr  for  maximum  gradient. 

2.2  Ti02  model  structures 

The  2-DPS  model  for  the  unfaceted  surface  used  a  (4x1)  supercell  with  six  layers,  giving 
48  Ti02  units  per  cell.  The  lateral  dimensions  were  chosen  to  give  a  low  coverage  of  adsorbed 
DMMP  (thus  avoiding  intermolecular  interaction)  while  still  being  computationally  affordable. 
The  surface  area  of  the  supercell,  containing  one  adsorbed  DMMP,  was  15.13  A  x  9.54  A  = 

144.2  A^,  and  the  coverage  was  1/8-monolayer  (one  DMMP  per  eight  surface  Tisc  sites).  A  slab 
thickness  of  six  layers  has  been  shown  [22,37]  to  be  adequate  for  convergence  of  the  surface  en¬ 
ergy  of  anatase  (100).  For  the  faceted  surface  the  supercell  was  increased  to  (4x2)  in  order  to  ac¬ 
commodate  the  facet.  After  removal  of  the  appropriate  atoms  (cf  Fig.  1),  the  supercell  comprised 
80  Ti02  units.  In  forming  the  facet,  equivalent  atoms  on  either  side  of  the  slab  were  removed  so 
as  to  maintain  inversion  symmetry.  Complete  relaxation  of  the  faceted  surface  will  be  seen 
(Sec.  3.1)  to  require  a  thicker  slab.  Thus  calculations  for  the  bare  faceted  surface  were  also  done, 
for  comparison,  for  a  thickness  of  twelve  layers  and  a  (1x2)  supercell.  Since  localized,  rather 
than  plane -wave,  basis  sets  are  used  one  can  study  an  isolated  2-DPS.  Hence  the  issue  of  the 
vacuum  gap  between  slabs  periodically  repeated  in  the  surface-normal  direction  does  not  arise. 

A  single  DMMP  molecule  was  placed  on  the  “top”  surface  of  the  2-DPS  supercell,  and  a 
second  was  placed  symmetrically  on  the  “bottom”  surface.  This  “two-sided”  approach  preserves 
the  center  of  symmetry  of  the  supercell,  thus  avoiding  the  production  of  a  dipole  potential,  and 
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has  been  shown  [43]  to  be  beneficial  in  2-DPS  studies  of  H2O  adsorption  on  rutile  (110).  The  ad¬ 
sorption  energy  per  DMMP  is  then  defined  as 

AEads  =  [E(Ti02+nDMMP)  -  E(Ti02)  -  nE(DMMP)  +  AE(BSSE)]/n 
where  n  =  2,  E(Ti02+nDMMP)  is  the  relaxed  total  DET  energy  with  adsorbed  DMMP,  E(Ti02) 
and  E(DMMP)  are  the  energies  of  bare  Ti02  and  free  DMMP  and  AE(BSSE)  is  the  counterpoise 
correction  for  basis  set  superposition  error.  Thus  AEads  <  0  indicates  an  exothermic  process.  All 
technical  parameters  (integration  accuracy,  etc.)  were  the  same  for  all  quantities  to  be  combined 
in  obtaining  AEads-  Specifically,  E(DMMP)  was  obtained  in  the  same  manner  as  the  other  energy 
terms;  namely,  by  optimizing  the  geometry  at  the  RHE/6-31G(d,p)  level  followed  by  a  single¬ 
point  B3EYP/6-31G(d,p)  calculation.  The  free  DMMP  molecule  was  placed  in  the  lowest-energy 
gas-phase  conformation  [44,45]  before  starting  the  optimization  process. 

3,  Results 

3.1  Relaxation  of  the  bare  surfaces 

Table  1  shows  the  lattice  parameters  obtained  in  the  present  RHP  calculation  for  bulk  ana- 
tase  and,  for  comparison,  bulk  rutile.  In  both  cases,  the  agreement  with  experiment  [46]  is  seen  to 
be  quite  good.  Table  2  shows  the  relaxation- induced  displacements  of  atoms  on  and  near  the  bare 
unfaceted  surface  in  comparison  to  previous  RHP  results  [37]  for  a  somewhat  thicker  slab,  toge¬ 
ther  with  corresponding  results  [22]  using  DPT  in  the  geometry  relaxation.  The  agreement 
among  the  three  sets  of  results  is  again  quite  good.  The  relaxation  is  essentially  identical  for  all 
equivalent  atoms  in  the  present  (4x1)  cell,  which  indicates  no  tendency  toward  reconstruction. 
Pigure  2  and  Table  3  show  corresponding  results  for  the  faceted  surface,  which  are  more  com¬ 
plex.  Por  the  2-DPS  described  above,  with  6  layers  before  faceting,  there  are  differences  in  the 
displacements  of  equivalent  atoms  having  different  positions  in  the  [001]  direction.  Por  example. 
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the  two  Layer  1  Tisc  sites  should  be  equivalent  (cf.  Fig.  lb)  but  exhibit  somewhat  different  dis- 
plaeements  during  relaxation.  However,  equivalent  atoms  at  different  positions  in  the  [100]  di- 
reetion  have  essentially  identieal  relaxations. 

These  anomalies  result  from  the  faet  that  a  thiekness  of  only  6  layers  is  insuffieient  for  eom- 
plete  relaxation  when  the  slab  is  thinned  by  faeeting.  This  is  shown  by  results  for  a  slab  with  12 
layers  before  faeeting  and  a  (1x2)  supereell.  Redueing  the  eell  dimension  in  the  [100]  direetion  is 
permissible  in  view  of  the  results  obtained  for  the  6-layer  slab.  For  the  thieker  slab  the  displaee- 
ments  of  equivalent  atoms  are  now  essentially  identieal,  but  these  are,  in  most  oases,  not  very 
different  from  those  for  the  thinner  slab.  The  signs  given  for  Ax  in  Table  3  are  signifioant.  Thus 
Oi(l)  and  0_l(2)  in  the  Layer  1  displaoe  away  from  the  faoet,  and  Oi(l)  and  Oi(2)  in  the  Layer  2 
displaoe  away  from  the  "ohannel".  In  Layer  6  of  the  12-layer  slab  (not  shown)  the  displaoements 
from  the  bulk  positions  are  all  <0.03  A  in  magnitude,  indioating  reasonably  good  oonvergence. 
Comparison  of  Figs.  1  and  2  shows  that  relaxation  leads  to  an  atomio-soale  roughening  of  the 
surface,  as  noted  previously  [22,37],  with  0||  (Tisc)  displacing  outward  (inward)  along  the  surface 
normal.  Finally,  the  density  of  states  for  either  surface  (not  shown)  revealed  no  surfaces  states  in 
the  band  gap,  in  agreement  with  previous  results  [10,17]  for  the  relaxed  unfaceted  surface. 

3.2  Molecular  adsorption  ofDMMP 
3.2.1  The  unfaceted  surface 

Figure  3  shows  the  relaxed  configuration  ofDMMP  molecularly  adsorbed  on  either  surface. 
For  the  unfaceted  surface  (Fig.  3a)  AEads  =  -21.0  kcal/mol  was  found  after  a  BSSE  correction  of 
+6.5  kcal/mol.  The  Ti— OP  and  TiO=P  bond  lengths  are  2.09  and  1.50  A  respectively,  and  the 
Ti — 0=P  bond  angle  is  138°.  There  appear  to  be  two  C-H — 02c  interactions  that  meet  the  H- 
bond  criteria  [47]  for  interatomic  distances  and  angles.  Eor  Bond  1  (2)  the  H — O  and  C — O  dis- 
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tances  are  2.28  and  3.29  A  (2.65  and  3.67  A)  respectively,  and  the  C-H — O  angle  in  either  case 
is  157°.  Other  possible  C-H— O  interactions  (based  on  the  H — O  distances)  involve  much  small¬ 
er  angles  generally  indicative  of  a  weaker  bond.  It  is  noted  that  AEads  for  the  (010)  surface  is 
close  to  2-DPS  results  obtained  elsewhere  [42]  for  molecular  adsorption  of  DMMP  on  the  ana- 
tase  (101)  and  rutile  (110)  surfaces  (-19.1  and  -17.6  kcal/mol,  respectively).  The  small  variation 
in  AEads  reflects  the  trend  in  surface  energies  given  above,  with  AEads  decreasing  in  magnitude 
with  increasing  surface  stability.  On  all  surfaces,  adsorption  reverses  the  inward  displacement  of 
the  Tisc  bonding  site  that  occurs  when  the  bare,  ideally-terminated  surface  relaxes.  This  can  be 
seen  clearly  in  Eig.  3a  as  an  outward  displacement  of  the  Tisc  adsorption  site  relative  to  the 
neighboring  Tisc  sites. 

Other  configurations,  differing  from  that  in  Eig.  3a  by  rotation  of  the  DMMP  by  90°  or  180° 
about  the  P=0  bond,  were  also  considered  as  starting  points  in  the  geometry  optimization.  These 
differ  as  to  which  combinations  of  CH3  groups  are  positioned  favorably  for  C-H— Oic  bond  for¬ 
mation.  At  the  RHF  level  used  in  geometry  optimization  the  90°  (180°)  rotations  relaxed  to  ener¬ 
gies  which  were  higher  by  about  2.3  (9.0)  kcal/mol  than  that  of  Eig.  3a.  Another  starting  geome¬ 
try,  differing  from  that  in  Eig.  3a  by  a  rotation  of  the  DMMP  by  180°  about  the  Ti — OP  bond, 
was  also  considered.  This  structure  makes  use  of  the  slightly  more  distant  O2C  sites  for  possible 
H-bonding.  However,  the  relaxed  structure  was  less  stable  by  ~15  kcal/mol  at  the  RHF  level. 

Still  other  rotations,  by  90°  about  the  Ti— OP  bond  or  by  180°  about  both  the  P=0  and  Ti— OP 
bonds  also  relaxed  to  structures  which  were  unstable  relative  to  that  in  Fig.  3a. 

The  fact  that  none  of  these  starting  structures  relaxed  to  the  apparently  most  stable  configur¬ 
ation  (Fig.  3a),  but  instead  remained  in  approximately  the  initial  conformation,  indicates  that 
there  is  a  significant  barrier  to  the  bond  rotations  needed  to  accomplish  the  transformation.  This 
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may  arise  from  the  need  to  break  one  or  more  C-H— 02c  bonds  in  order  to  aecomplish  the  rear¬ 
rangement.  The  existenee  of  this  barrier  made  it  neeessary  to  eonsider  several  different  initial 
struetures  in  the  geometry  optimization. 

Finally,  dative  bonding  between  a  Tisc  and  a  methoxy  O  atom  is  also  a  possible  mode  of  ad¬ 
sorption.  This  has  been  eonsidered  elsewhere  [42]  for  the  rutile  (110)  surfaee  and  been  found  to 
give  a  negligible  AEads-  In  the  present  ease,  an  attempt  was  made  to  adsorb  DMMP  on  the  unfa- 
eeted  surfaee  via  a  Tisc — O  bond  to  a  methoxy  O  atom.  At  the  RHF  level  used  in  geometry  opti¬ 
mization  the  resulting  strueture  was  less  stable  than  that  in  Fig.  3a  by  '-'25  keal/mol,  indicating  an 
unstable  configuration.  Likewise,  an  initial  structure  with  a  Tisc— 0=P  bond  together  with  a 
Tisc — O  bond  to  a  methoxy  at  an  adjacent  Tisc  site  relaxed  to  one  with  only  the  Tisc — 0=P  bond. 
3.2.2  The  faceted  surface 

For  the  faceted  surface  (Fig.  3b),  DMMP  was  placed  symmetrically  on  both  slab  faces  at  the 
Ti5c(l)  site  of  Layer  2  shown  in  Lig.  2.  The  slab  at  this  point  is  four  layers  thick,  which  should 
suffice  for  a  nearly-converged  ALads,  and  is  preferable  to  the  Layer  2  Ti5c(2)  site  where  the  slab  is 
five  layers  thick.  Previous  work  on  the  adsorption  of  CO  [48]  or  H2O  [49,50]  on  rutile  (110) 
shows  that  a  thin  2-DPS  with  an  even  (odd)  number  of  Ti  layers  tends  to  slightly  underestimate 
(significantly  overestimate)  the  fully-converged  magnitude  of  ALads-  Similar  odd/even  oscilla¬ 
tions  in  the  (100)  surface  energy  are  seen  for  thin  anatase  2-DPS  slabs  [37]. 

Lor  the  faceted  surface,  ALads  =  -21.3  keal/mol  was  found  using  the  same  AL(BSSL)  as 
above.  The  Ti — OP  and  TiO=P  bond  lengths  are  2.04  and  1.49  A  respectively,  and  the  Ti— 0=P 
angle  is  140°.  All  results  are  close  to  those  given  above  for  adsorption  on  the  unfaceted  surface. 
The  dashed  line  indicates  one  of  several  possible  H-bonding  interactions,  based  on  the  H — O  dis¬ 
tances  which  are  in  the  range  of  2.39  to  3.00  A.  The  indicated  bond  is  believed  to  be  the  strong- 
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est,  with  O — H  and  C — O  distances  of  2.39  and  3.29  A  respectively  and  a  C-H — O  angle  of 
140°.  Other  interactions  have  smaller  C-H — O  angles  and/or  longer  H— O  distances  and  are 
therefore  probably  weaker  [47].  It  is  possible,  however,  that  the  cumulative  effect  of  several 
weak  H-bonds  might  make  a  significant  contribution  to  the  total  AEads-  A  quantitative  discussion 
of  H-bonding  requires  the  characterization  of  bond  critical  points  using  the  Atoms  in  Molecules 
theory  [51]  which  is  not  available  at  present. 

Alternative  starting  structures,  differing  from  that  in  Fig.  3b  by  rotations  of  70  or  140“  about 
the  P=0  bond,  were  also  considered  in  order  to  position  various  CH3  groups  for  possible 
H-bonding  to  O2C  sites  in  Layers  1  and  2.  However,  these  were  less  stable  by  about  3.5  and 
7.2  kcal/mol  respectively  at  the  RHF  level  used  in  geometry  optimization.  As  in  the  case  of  the 
unfaceted  surface,  neither  of  these  initial  structures  relaxed  to  the  most  stable  configuration  dur¬ 
ing  geometry  optimization  which  suggests  a  significant  barrier  to  the  necessary  bond  rotation. 
Large  rotations  about  the  Ti — OP  bond  were  also  considered  but  were  hindered  by  steric  interac¬ 
tions  between  the  CH3O  groups  and  the  "step  edge"  between  Layers  1  and  2. 

As  noted  above,  both  The  sites  in  Layer  2  (Fig.  2)  are  equivalent  on  the  semi-infinite  surface 
but  can  differ  for  a  thin  2-DPS.  Adsorption  at  the  Ti5c(2)  site  of  Layer  2  was  evaluated  for 
DMMP  in  the  same  configuration  as  at  the  Ti5c(l)  shown  in  Fig  3b.  The  bond  lengths  and  angles 
in  the  relaxed  structure,  including  those  of  the  putative  dominant  H-bond,  were  virtually  identical 
to  those  found  for  adsorption  at  Ti5c(l),  but  AEads  was  larger  in  magnitude  by  ~4.3  kcal/mol  at 
the  RHF  level  used  in  geometry  optimization.  On  the  basis  of  the  above  discussion,  this  is  as¬ 
sumed  to  represent  an  overestimate.  Adsorption  outside  the  facet  (i.e.,  at  a  Layer  1  Tisc  site  in 
Fig.  3b)  was  compared  with  adsorption  inside  the  facet  (i.e.,  at  the  Layer  2  Tisc  site  as  shown  in 
Fig.  3b).  The  starting  geometry  was  the  same  as  those  identified  above  as  the  most  favorable  and 
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shown  in  Figs.  3a, b.  At  the  RHF  level  the  resulting  AEads  was  3.9  Keal/mol  smaller  in  magnitude 
than  for  adsorption  inside  the  faeet.  The  overall  eonclusion  is  that  the  results  do  not  depend 
greatly  on  the  ehoiee  of  adsorption  site  but  that  the  strueture  in  Fig.  3b  is  the  most  favorable  for 
moleeular  adsorption  on  the  faeeted  surfaee. 

3.3  Dissociative  adsorption  of  DMMP 

Figure  4  shows  models  for  a  possible  first  step  in  the  dissociation  of  DMMP  on  surfaces 
with  and  without  faceting.  The  structure  is  obtained  by  breaking  a  PO-CH3  bond,  moving  the 
CH3  to  an  O2C  site  and  then  forming  a  (CH3)(CH30)P(-0-)2  bridge  between  two  Tisc  sites.  This  is 
only  a  formal  description  and  is  not  intended  to  represent  the  actual  reaction  mechanism.  The  end 
product  is  consistent  with  the  results  of  x-ray  photoemission  spectroscopy  (XPS)  and  tempera¬ 
ture-programmed  desorption  (TPD)  experiments  for  DMMP  adsorbed  at  room  temperature  on 
OH-free  single-crystal  rutile  (110)  surfaces  prepared  and  maintained  in  UHV  [29,30].  These  data 
show  mainly  molecular  adsorption.  The  relatively  small  amount  of  decomposition  that  does  oc¬ 
cur  leads  to  the  desorption  of  CH3,  but  not  CH3O,  in  TPD.  The  yield  of  CH3  is  distinct  from  the 
cracking  of  desorbed  DMMP  itself  and  is  consistent  with  CH3  adsorbed  at  a  bridging  O2C  site  as 
opposed  to  a  terminal  CHs-O-Ti  site.  Desorption  of  CH3O  from  such  a  terminal  site  would  re¬ 
form  a  Tisc;  whereas,  desorption  of  CH3  would  yield  an  unstable,  doubly-unsaturated  O  atom  un¬ 
less  a  Tisc  site  were  available  for  forming  an  O2C  bridge.  However,  these  data  were  obtained  for 
rutile  (110)  and  might  not  necessarily  apply  directly  to  anatase  (010). 

A  thorough  examination  of  all  a  priori  possible  dissociation  products  was  not  performed 
since  the  goal  here  is  to  determine  the  extent  to  which  faceting  affects  dissociation,  and  the  rela¬ 
tive  energies  for  molecular  vs.  dissociative  adsorption  for  a  plausible  structure.  Several  other  dis¬ 
sociation  structures  are  considered  elsewhere  [42]  for  DMMP  on  rutile  (110)  and  anatase  (101). 
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The  starting  configuration  was  based  on  the  most-favorable  structures  for  molecular  adsorption 
discussed  above  and  shown  in  Fig.  3.  For  the  unfaceted  surface,  one  alternative  structure  was 
considered  in  which  the  CH3  was  placed  on  one  of  the  two  equivalent  O2C  sites  which  are  closest 
to  the  bridge  (i.e.,  one  step  to  the  right  of  that  occupied  in  Fig.  4b).  This  alternative  site  is  marked 
by  an  arrow  in  Fig.  4b.  The  OCH3  on  the  bridge  was  also  rotated  90°  about  the  P-OCH3  bond  so 
as  to  move  the  CH3  away  from  the  O2C-CH3,  thus  minimizing  steric  repulsion  between  the  two 
CH3  groups  and  allowing  C-H — O  bonding  between  the  O2C-CH3  and  the  O  atom  in  the  P-O-CH3 
group.  At  the  RHF  level  the  geometry-optimized  AEads  was  negligibly  different  from 
(~1  kcal/mol  larger  than)  that  for  the  structure  in  Fig.  4b. 

The  BSSE-corrected  AEads  is  -32.7  and  -22.7  kcal/mol,  respectively,  for  the  unfaceted  and 
for  the  faceted  surface.  These  are  larger  than  the  energies  given  above  for  molecular  adsorption 
but  only  slightly  so  for  the  faceted  surface.  This  indicates  that  dissociation  is  thermodynamically 
favored  to  a  degree  which  depends  on  faceting.  As  noted  above,  faceting  exposes  the  more  stable 
(Oil)  surface,  which  is  isostructural  with  the  (101).  Consistent  with  this,  a  2-DPS  calculation 
done  elsewhere  [42]  for  the  dissociative  adsorption  of  DMMP  on  the  anatase  (101)  surface,  in  a 
structure  like  that  studied  here,  gave  a  BSSE-corrected  AEads  of  -25.0  kcal/mol  which  is  close  to 
the  present  result  for  the  faceted  (010).  Although  no  relevant  data  are  available  for  completely 
OH-free  anatase,  DMMP  adsorption  on  OH-free  rutile  (1 10)  at  room  temperature  is  largely  mo¬ 
lecular  [29,30].  If  the  same  is  true  for  anatase  (010)  then  the  present  results  suggest  that  dissocia¬ 
tion  under  these  conditions  is  kinetically  limited.  This  is  consistent  with  experimental  results  for 
partially-hydroxylated  Ti02  powder  [31,32]  which  indicate  that  surface  OH  is  directly  involved 
in  the  decomposition  of  molecular  DMMP,  possibly  via  a  nucleophilic  attack  on  the  P  atom  by 
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the  O  atom  of  a  Ti-OH  group,  and  that  this  can  occur  below  room  temperature.  In  the  present  in¬ 
terpretation  the  Ti-OH  is  necessary  to  overcome  the  kinetic  barrier  to  decomposition. 

4,  Conclusions 

This  is,  to  our  knowledge,  the  first  computational  study  of  the  faceted  anatase  (010)  (or, 
equivalently,  (100))  surface  either  with  or  without  an  adsorbate.  On  the  bare  surfaces,  similar  re¬ 
laxations  occur  for  corresponding  atoms  with  and  without  faceting.  These  consist  mainly  of  in¬ 
ward  (outward)  displacements  of  Tisc  (Tiec)  atoms  and  outward  displacements  of  Osc  atoms.  Sur¬ 
face  02c  atoms  displace  mainly  parallel  to  the  surface.  The  energies  and  structures  for  molecular 
adsorption  of  DMMP  are  similar  in  either  case  and  are  also  similar  to  those  [42]  for  rutile  (110) 
and  anatase  (101)  surfaces.  Adsorption  occurs  preferentially  via  a  Tisc— 0=P  dative  bond  with 
evidence  of  C-H — O  bonding  between  methoxy  groups  and  02c  atoms.  The  results  indicate  that 
molecular  adsorption  of  DMMP  on  HSA  Ti02  powders,  for  which  the  anatase  (010)  surface  is  a 
major  constituent,  can  be  understood  without  necessarily  taking  account  explicitly  of  the  facet¬ 
ing.  As  a  corollary,  the  results  also  suggest  that  data  for  molecular  adsorption  of  DMMP  on  the 
easily-obtained  rutile  (110)  surface  should  be  transferable  to  the  less  common  anatase  (010)  surf¬ 
ace.  However,  this  applies  only  to  the  pristine  surfaces  considered  here  which  are  free  of  defects 
and  adsorbed  OH.  The  situation  may  change  when  such  imperfections  are  introduced. 

A  possible  first  step  in  the  dissociative  adsorption  of  DMMP  on  OH-free  surfaces  was  de¬ 
fined  on  the  basis  of  experimental  data  [29,30]  for  OH-free  rutile  (110)  single-crystal  surfaces. 
On  the  non-defective  and  OH-free  surface  the  product  consists  of  an  O2C-CH3  bond  at  a  Ti-02c-Ti 
site  and  a  (CH3)(CH30)P(-0-)2  bridge  between  two  Tisc  sites.  AEads  for  this  process  indicates  that 
it  is  thermodynamically  favored  over  molecular  adsorption  but  to  a  degree  that  depends  on  facet¬ 
ing.  Faceting  exposes  the  more  stable  (011)  (or,  equivalently,  (101))  surface  and  reduces  the 


15 


magnitude  of  AEads  for  dissociative  adsorption,  shifting  it  toward  the  result  obtained  for  the 
(Oil).  Experiments  on  fully  OH-free  anatase  (010)  surfaces  are  needed  in  order  to  determine 
whether  dissociation  is  kinetically  accessible  without  the  intervention  of  Ti-OH  sites. 
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Table  1 


Observed  and  ealculated  Ti02  lattice  parameters^ 


rutile*’ 

anatase‘^ 

expt. 

calc. 

expt. 

calc. 

a 

4.5937 

4.564 

3.7845 

3.781 

c 

2.9587 

2.987 

9.5143 

9.536 

u 

0.30478 

0.3060 

0.20806 

0.2063 

Volume 

62.435 

62.224 

136.268 

136.355 

^  Experimental  values  are  from  neutron  diffraction  data,  ref  46.  The  lattice  constants  (a,  c)  are  in 
Angstroms,  and  u  is  a  fraction  of  a  lattice  constant.  The  volume  of  the  crystallographic  unit  cell 
(which  comprises  two  primitive  cells  in  the  case  of  anatase)  is  given  in  A^. 

’’  The  rutile  Ti  positions  are  at  (0,0,0)  and  (14, 14, 14).  The  O  positions  are  at  ±(u,  u,  0), 

(-U+14,  u+14, 14)  and  (u+14,  -u+14, 14). 

The  anatase  Ti  positions  are  at  ((0,  0,  0),  (14, 14, 14))  +  ((0,  0,  0),  (0, 14,  Va)).  The  O  positions  are 
at  ((0,  0,  0),  (1/2,  /2,  /))  +  ((0,  0,  u),  (0,  /2,  U+/4),  (1/2,  0,  -m+Ya),  (/,  /,  -u+Z)) 
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Table  2 


Displacements,  relative  to  the  ideally-terminated  surfaee,  resulting  from  relaxation  of  the  unfa¬ 
ceted  anatase  (010)  surfaee^^ 


atom'’ 

this  work‘’ 

ref  37^ 

ref  22® 

Tisc 

0.02; -0.15 

0.04; -0.14 

0.02; -0.16 

Layer  1 

0| 

0.02;  +0.18 

0.04;  +0.18 

0.04;  +0.18 

Ox 

0.17;  0.00 

0.16;  -0.01 

0.16;  +0.02 

Ti6c 

0.00;  +0.12 

0.00;  +0.11 

0.01; +0.17 

Layer  2 

0| 

0.00;  +0.07 

0.00;  +0.07 

0.01; +0.10 

Ox 

0.12; -0.04 

0.11; -0.03 

0.13; -0.03 

^  All  displaeements  are  in  Angstroms.  Values  are  given  as  |Ax|;  Az  where  Ax  is  the  displaeement 
in  the  [001]  direetion,  and  Az  is  the  displaeement  along  the  surfaee  normal.  A  positive  (negative) 
Az  is  outward  (inward).  Layer  1  is  the  surfaee  layer.  For  Oi,  Ax  is  in  a  direetion  away  from  the 
"ehannel"  (ef.  Fig.  1).  The  anatase  (100)  and  (010)  surfaces  are  isostructural. 

’’  0||  refers  to  the  O  atom  bonded  to  three  Ti  atoms  in  the  same  plane  whieh,  at  the  surfaee,  is  an 
Osc  site.  0_L  refers  to  the  O  atom  with  bonds  to  Ti  in  different  planes  whieh,  at  the  surfaee,  is  an 
02c  site.  See  Fig.  1. 

Unfaeeted  (4x1)  2-DPS  model  with  six  layers  after  geometry  relaxation  as  deseribed  in  the  text. 
(lOO)-(lxl)  2-DPS  model  with  eight  layers  relaxed  in  an  RHF  ealeulation  with  ECPs  for  Ti  and 
O  (see  ref  37). 

®  (lOO)-(lxl)  2-DPS  with  six  Ti  layers  relaxed  in  a  DFT  ealeulation  using  Car-Parrinello  molecu¬ 
lar  dynamics  (see  ref  22). 
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Table  3 


Displacements,  relative  to  the  ideally-terminated  surfaee,  resulting  from  relaxation  of  the  faeeted 
anatase  (010)  surface^^ 


Layer  1 

6  Layers  12  Layers 

Layer  2 

6  Layers  12  Layers 

atom’’ 

Ax;  Az 

Ax;  Az 

atom’’ 

Ax;  Az 

Ax;  Az 

Ti5e(l) 

-0.05;  -0.08 

-0.02; -0.10 

Ti5c(l) 

+0.06;  -0.25 

0.00;  -0.16 

Tisc  (2) 

-0.01; -0.13 

+0.02; -0.10 

Tisc  (2) 

0.00;  -0.11 

0.00;  -0.16 

0||(1) 

-0.05;  +0.25 

-0.02;  +0.23 

Ti6c  (1) 

+0.04;  +0.08 

+0.02;  +0.15 

011  (2) 

0.00;  +0.21 

+0.03;  +0.23 

Ti6c  (2) 

-0.04;  +0.17 

-0.02;  +0.15 

Ou(l) 

+0.13;  0.00 

+0.16;  +0.04 

0||(1) 

+0.05;  +0.09 

+0.01; +0.20 

Ou(2) 

-0.19;  +0.05 

-0.16;  +0.04 

011  (2) 

-0.02;  +0.25 

-0.01; +0.20 

0  |  (3) 

+0.08;  0.00 

+0.04;  +0.06 

011  (4) 

-0.05;  +0.07 

-0.04;  +0.06 

Ou(l) 

+0.24;  +0.01 

+0.19;  -0.02 

Ou(2) 

-0.19;  0.01 

-0.18; -0.01 

Ou(3) 

+0.13; +0.01 

+0.12;  +0.01 

Ou(4) 

-0.14; -0.01 

-0.12;  +0.01 

^  See  Table  2,  footnote  a.  "6  Layers"  refers  to  the  2-DPS  with  a  (4x2)  supereell  and  6  layers  (be¬ 
fore  introduction  of  the  faeet).  "12  Layers"  refers  to  the  2-DPS  with  a  (1x2)  supercell  and  12  lay¬ 
ers  (before  faeeting).  See  See.  3.1  for  further  details. 

’’  See  Table  2,  footnote  b,  and  Fig.  2. 
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Figure  captions 

Fig.l  The  anatase  (010)  surface  with  the  bulk  unit  cell  shown  by  the  dashed  lines.  Blue  (red) 
spheres  are  Ti  (O).  (a)  Shows  the  ideally-terminated  surface,  and  (b)  shows  a  (1x2)  faceted  struc¬ 
ture  (ref.  23).  Both  structures  are  shown  before  relaxation.  The  different  models  are  not  to  scale. 
The  encircled  section  in  (a)  shows  what  is  removed  in  forming  the  facet.  Various  surface  sites 
are  labeled  (see  text),  and  the  region  referred  to  as  the  "channel"  in  the  text  is  indicated.  The 
plane  shown  in  (b)  is  the  (Oil),  which  is  isostructural  with  the  (101). 

Fig.  2  The  faceted  surface  after  relaxation.  In  correspondence  to  Table  2,  the  brackets  indicate 
the  layer  numbers,  and  the  individual  atoms  are  labeled  by  type.  For  clarity,  Ti  (O)  atoms  are  la¬ 
beled  on  the  "top"  ("bottom")  of  the  slab.  The  Layer  1  Ti5c(2)  is  obscured  by  an  O  atom.  The 
Layer  2  Tisc  site  shown  italicized  in  blue  is  the  one  used  for  DMMP  adsorption  within  the  facet 
(see  text).  The  [001]  and  [010]  vectors  define  the  positive  directions  for  Ax  and  Az,  respectively, 
for  the  "top"  surface  atoms  (cf  Table  2  and  discussion  in  text).  The  [001]  and  [010]  directions 
are  reversed  on  the  “bottom”  surface 

Fig.  3  Sections  of  the  (a)  unfaceted  and  (b)  faceted  surface  after  relaxation  with  adsorbed 
DMMP  in  the  lowest-energy  conformations.  For  clarity,  only  sections  of  the  entire  slab  super¬ 
cells  are  shown.  The  dashed  lines  show  H-bonding  interactions  identified  as  such  on  the  basis  of 
the  02c — H  and  Oic — C  distances  and  the  02c — H-C  angle  (see  text).  In  (b)  the  adsorption  site  is 
the  one  designated  Layer  2  Ti5c(l)  in  Fig.  2.  Note  that  the  orientation  in  (b)  is  reversed  relative  to 
that  in  Fig.  2.  The  green,  black  and  gray  spheres  represent  P,  C  and  H  respectively. 

Fig.  4  Similar  to  Fig.  3  but  showing  structures  considered  for  the  first  step  in  DMMP  dissocia¬ 
tion.  The  upper  panel  shows  the  unfaceted  surface  viewed  along  the  (a)  [100]  and  (b)  [001]  di¬ 
rection.  (c)  Shows  the  faceted  surface  viewed  along  the  [100]  direction.  Possible  H-bonding  in- 
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teractions  have  been  omitted  for  clarity.  The  arrow  in  Fig.  4b  marks  an  alternative  location  for 
the  CH3  fragment  which  is  discussed  in  Sec.  3.3. 
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